Kinetics of bond formation in cross-linked gelatin gels 
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In chemical cross-linking of gelatin solutions, two different time scales affect the kinetics of the gel 
formation in the experiments. We complement the experimental study with Monte Carlo numerical 
simulations of a lattice model. This approach shows that the two characteristic time scales are related 
to the formation of single bonds cross linker-chain and of bridges between chains. In particular their 
ratio turns out to control the kinetics of the gel formation. We discuss the effect of the concentration 
of chains. Finally our results suggest that, by varying the probability of forming bridges as an 
independent parameter, one can finely tune the kinetics of the gelation via the ratio of the two 
characteristic times. 

PACS numbers: 



I. INTRODUCTION 



Among biopolymers, gelatin gels have received great 
attention [l[ because of their numerous applications 
in pharmaceutical, photographic and food industries. 
When a semi-diluted gelatin solution is cooled below 
room temperature, the coils start to form triple helices 
and progressively a connected network is built. The triple 
helices are reminiscent of the structure of native collagen, 
which gave origin to gelatin by means of a denaturation 
process. The gel is thermoreversible, i.e. by raising the 
temperature the sol state is recovered. Biodiversity due 
to chemical composition of the native collagen, molecular 
weight distribution, solution properties such as concen- 
tration or pH, influences the temperature of helix for- 
mation in the physical gel Q. The shear modulus shows 
universal scaling behavior with a critical exponent / close 
to 2 versus the distance from the critical concentration 
of helices Q. On the other hand, gelatin solutions show 
an even richer phenomenology since chemical gelation 
or a combination of chemical and physical gelation can 
be observed. In fact, if the system is kept above the 
helix formation temperature, amino-acids present along 
the gelatin chain can react with cross-linking molecules 
added to the solution. In this case helices cannot form 
and a permanent network appears due to cross-links be- 
tween reactant and chains, leading to the onset of an 
elastic response. Recently, extended studies have been 
performed on systems of gelatin in solution with bisvinyl- 
sulphonemethyl (BVSM) reactant [H, able to establish 
double covalent bonds with the lysine, the hydroxy ly- 
sine and possibly with other amine groups of gelatin 
chains. The chemical reaction is schematically shown in 
Figffl The influence of various parameters, as gelatin or 
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FIG. 1: Reaction leading to bond formation between an amine 
group and BVSM 



reagent concentration and solution pH, on bond forma- 
tion was investigated. For instance, increasing the pH ac- 
tivates more amine groups able to react with BSVM along 
the gelatin chains. Microcalorimetry measurements were 
able to monitor the development of the chemical reaction 
in time by detecting the exothermic enthalpy change dur- 
ing the formation of C — N bonds. Then the kinetics of 
cross-link formation was found to follow a double expo- 
nential decay with two characteristic times, whereas a 
simple exponential decay was detected at low pH. 

However, when counting the number of cross-links 
binding BVSM and gelatin, the method could not dis- 
criminate between bonds established by free reactants 
with a chain and bonds leading to a bridge between two 
gelatin chains, nor else bonds leading to a loop within a 
chain. This lack of information on the kinetics leading to 
the gel formation crucially affects the characterization of 
the gel structure and therefore its mechanical properties. 
Finally it may reflect onto the location of the gelation 
threshold and the determination of the critical exponent 
of the elastic modulus (the critical behavior of the shear 
modulus was measured at low frequency, giving a critical 
exponent / = 3.4 ± 0.3 [4|, close to the expected value for 
the vulcanization of long chains). As a consequence, a 
deeper comprehension of the bond formation kinetics is 
essential, requiring alternative investigations. In partic- 
ular the primary question is to understand how the two 
time scales controlling the kinetics depend on the forma- 
tion of single-bonds and bridges between the cross-linkers 
and the chains, or else to loops within the chains. More- 
over it would be crucial to understand how these time 
scales are related to the properties of the gelatin solution 
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(concentration, pH...) and of the cross-linking molecules 
(concentration, reactivity...). 

In this paper we analyze the role of concentration and 
reactivity of cross-linking molecules in the kinetics of 
bond formation, by complementing the experimental ob- 
servations with a numerical study. In gelling systems 
numerical approaches to the study of rheological and dy- 
namical properties have revealed to be extremely useful 
for a better understanding of experimental data [!, H, 0] • 
Both Monte Carlo and molecular dynamics simulations 
have been applied in the last years to the study of dif- 
ferent chemical gelation processes 0, [8, 9] . Here we have 
used Monte Carlo simulations on the cubic lattice of 
a simple model to analyze the kinetics of bond forma- 
tion in chemical cross-linking of a gelatin solution. We 
have considered a solution of polymer chains at differ- 
ent concentrations. Reactant monomers can diffuse in 
the solution forming bonds with the active sites along 
the chains, producing the cross-linking. Within this ap- 
proach we have followed the kinetics of the gel formation 
varying the gelatin concentration, the cross-linker con- 
centration and its bonding probability (i.e. reactivity). 
Our data reproduce extremely well the experimental find- 
ings. They show that the two time scales detected in 
the experiments correspond respectively to the average 
time of forming single bonds reactant-chains and bridges 
chains-chains via cross-linkers. The ratio of these two 
characteristic times controls the kinetics of the bond for- 
mation: Varying the concentration and the cross-linker 
reactivity strongly affect this ratio and therefore the ki- 
netics of the gelation process. 

The paper is organized as follows: In sect[TI] the re- 
sults of microcalorimetry measurements are presented, 
whereas in section Mil the numerical study is described 
and the results on bond formation are discussed. Finally 
in section IIVI the kinetics of bond formation is analyzed 
and the concluding remarks are given in section fVl 



II. MICROCALORIMETRY MEASUREMENTS 

In experiments, the kinetics of the reaction between 
amine groups of gelatin chains and BVSM has been mon- 
itored by microcalorimetry measurements. The gelatin 
sample is a photographic grade of gelatin extracted from 
lime processed ossein with an average molecular weight 
M w ~ lQbZOOg/mole, an index of polydispersity I p — 
2.06 and an isoelectric point pi = 4.9. The granules con- 
tain approximately 10% humidity and the concentrations 
are corrected accordingly. The BVSM can create cova- 
lent C — N bonds with amine groups of gelatin chains, 
so that a permanent network is formed at T > 40°C, 
where the triple helices structure of gelatin gels does not 
form. D. Hcllio-Scrughetti and M. Djabourov [|[ ana- 
lyzed the relation between elastic properties and system 
parameters was analyzed. The exothermal reaction be- 
tween amine groups and BVSM has been monitored by 
measuring the enthalpy change. We have performed sev- 
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FIG. 2: Released heat during the chemical reaction be- 
tween gelatin and reactant, for a solution with C ge i = 12%, 
Cbvsm = 0.15%, pH = 6.7 and T = 40°C 

eral experiments at temperature T = 40°C, using so- 
lutions with different concentrations of gelatin and re- 
actant, for different values of pH. In FigJ5] we plot the 
released heat Q(t) as a function of time for a solution 
with gelatin concentration C ge i — 12%g/cm 3 , BVSM 
concentration Cbvsm — 0.15%<?/cm 3 and pH = 6.7. 
Normalizing the released heat by the enthalpy change 
AH = —AQkJ/mol due to the formation of one C — N 
bond, the curve represents at any time the total number 
of bonds formed between gelatin chains and BVSM. By 
writing Q(t) — A(l — f(t)), where A is a dimensional co- 
efficient proportional to AH, we introduce the function 
f(t) which represents the fraction of bonds that remain 
to form at time t. In FigJ3] f(t) is plotted as a func- 
tion of time. Data have been fitted with the sum of two 
exponentials: 

f(t) = A x exp(-t/rD + A 2 exp(-t/r 2 m ) (1) 

with r|" = 520s, r 2 m = 9000s so that r^/rf = 17.31 
(the apex "m" is an abbreviation for "microcalorime- 
try"). It is worth to notice that microcalorimetry experi- 
ments do not allow to discriminate between single-bonds 
and bridges or loops within a chain, where only bridges 
between different chains contribute to the increase of con- 
nectivity in the system. For this reason computer simula- 
tions are a fundamental step for a deeper understanding 
of the kinetics of bonds formation. 



III. MODEL AND NUMERICAL STUDY 

We have performed Monte Carlo simulations on a cu- 
bic lattice of a system made of bi- functional monomers, 
i.e. the reactant and linear chains, which are represented 
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FIG. 3: The fraction of bonds of the BVSM that remain to 
form as a function of time, for a solution with C ge i — 12%, 
Cbvsm = 0.15, pH = 6.7 and T = 40°C. The continuous 
line is the fit with Eq.Q. 



by a sequence of n = 10 linked monomers. One monomer 
of the chain models a Kuhn segment [lj| , and therefore 
represents more units. The length of a Kuhn segment in 
a gelatin chain has been measured to be of the or- 
der of 40 A, corresponding to about 10 amino-acids. As 
compared to the experiments, our chains correspond to 
shorter gelatin chains, containing only about 100 amino- 
acids. Each monomer occupies simultaneously the eight 
sites of the lattice elementary cell and, to take into ac- 
count the excluded volume interaction, two occupied cells 
cannot have any site in common. Some monomers along 
the chain are active sites which may bind to the reactant 
in order to compose complex clusters of chains leading 
to the formation of a gel. The active sites are tetra- 
functional. Two bonds are formed with the neighbors 
along the chain and two are not saturated at the begin- 
ning of the simulation. The number of active sites per 
chain, n asi corresponds to a fixed pH of the solution. In 
fact, in experiments the increasing of the pH activates 
more amine groups able to react with the BVSM along 
the chain, therefore in simulations n as could be varied to 
study the effect of the pH. Although the number of amine 
groups in a gelatin chain actually linked to reactant can- 
not be measured experimentally, it is estimated that at 
most a fraction of 20% can react. Therefore we have 
performed most of the simulations for n as = 5, which 
corresponds to a fraction of 10% of active amino-acids in 
our chain. 

Chains are randomly distributed on the lattice and dif- 
fuse via random local movements. The excluded volume 
interaction and the self avoiding walk condition for poly- 
mer clusters restrict the possibility of monomer move- 
ments: to satisfy these two requirements the bond lengths 
vary into a set of permitted values according to bond- 



fluctuation dynamics [12| . On a cubic lattice the allowed 
bond lengths are 1 — 2, \/5, y/6, 3, VTO. At each Monte 
Carlo step the time is increased by St = 1 and one ran- 
dom move is selected on average for each monomer: if 
the move satisfies the bond-fluctuation dynamics and ex- 
cluded volume conditions, it is executed, otherwise it is 
rejected. These simple laws for local movements give rise 
to a dynamics which takes into account the main features 
of the real dynamics of polymer molecules [12] • 

After chains have diffused and reached equilibrium, we 
add the reactant to the system and let the solution diffuse 
towards the stationary state. Due to the diffusion of 
cross-linkers and chains, when a reactant finds a nearest 
neighbor unsaturated active site, a bond may form. The 
process goes on until all the possible bonds are formed. 

The bond formation may request to overcome a free 
energy barrier depending on the nature of the so- 
lution, the active sites and the reactant. In particular, 
it may depend on some specific local orientation of the 
molecules, some restriction on the value of the angle be- 
tween two bonds, due to the rigidity of the C — N link, or 
else may be affected by variations of the effective reactiv- 
ity of the cross-linker. In our model we have taken into 
account these effects in the following way: The first bond 
of a reactant monomer is formed along lattice directions 
as soon as there is a neighboring active site. The second 
bond is formed with probability pb < 1, since a reactant 
monomer is expected to have less chance to react when is 
already bonded to a chain, compared to when it is free. 
In the same spirit of reaction limited aggregation models 
[3] , Pb is an independent parameter which influences the 
time of formation of bridges between gelatin active sites. 
The value of pb should be determined by the features 
of cross-linking reagent. Moreover, although varying the 
bridge probability pb does not affect the gelation transi- 
tion, it has a crucial effect on the velocity of the reaction, 
which can be easily enlightened in the numerical simula- 
tion as discussed later. 

We have performed numerical simulations of the model 
for different lattice sizes (L = 50,100,200), where the 
unit length is the lattice spacing a = 1, with periodic 
boundary conditions. The chain concentration C and 
the cross-linker concentration C r are defined as the ra- 
tio between the number of monomers/reactant and the 
maximum number of monomers N max = L 3 /8 in the sys- 
tem. Using the percolation approach we identify the gel 
phase as the state in which there is a percolating clus- 
ter, which spans the whole system [l5| . For a fixed set 
of parameters we generate a number of configurations of 
the system and monitor the reaction. In order to locate 
the gelation transition we analyze the percolation proba- 
bility II, defined as the fraction of configurations leading 
to a percolating cluster, and we identify the transition 
with the line II = 0.5 [16|. We have determined a qual- 
itative phase diagram (Fig|4]) by varying the chain and 
cross-linker concentrations, C and C r respectively, for 
a fixed n as — 5. In experiments the total amount of 
reactant has been consumed at the end of the reaction 
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FIG. 4: (Color online) The phase diagram, obtained by plot- 
ting the percolation probability II as a function of chain and 
reactant concentration C and CV, respectively, using a color 
scale, reported in legend. The spanning probability has been 
averaged over 30 independent configurations of a sample of 
size L = 100. The number of active sites per chain n a3 = 5 
is kept constant and p& = 1. The percolation line can be 
identified with the locus II = 0.5. 

process, i.e. the amount of reactant is much lower than 
the amount of active sites. The reaction stops when all 
the reactant are linked to amine groups and the system 
is in the gel phase. Moreover the experimental system 
[3| is investigated at gelatin concentrations C ge i above 
the overlap concentration C* el = 0.005 g/cm~ 3 . In order 
to reproduce the experimental conditions the crosslinker 
concentration has been fixed at G = 0.025, which cor- 
responds to C r < C at the sol-gel transition. Under 
this condition, in simulations the gel phase is located at 
concentrations C above the overlap concentration which 
in our system is C* ~ 0.017. This choice of parameters 
guarantees that at the end of the reaction the system is 
in the gel phase. 

With finite size scaling analysis we have obtained the 
percolation threshold, that for the case C r — 0.025 is 
C c = 0.10 ± 0.05, the critical exponents v = 0.9 ± 0.1 for 
the percolation connectivity length £ (~ \C — C c \~ v ) and 
7 = 1.78±0.10 for the mean cluster size \ (~ \C— C c | -7 ). 
These results are in good a gree ment with the random per- 
colation critical exponents [16j. The random percolating 
cluster is characterized by a fractal structure: its mass 
M, i.e. the number of monomers, scales with its linear di- 
mension £ with a power law behavior £ D , where D ~ 2.5 
[l6[ . The structure of the formed network may depend on 
model parameters and influences the rheological response 
of the system. 



IV. KINETICS OF BOND FORMATION 

In simulations we have investigated the behavior of the 
number of bonds formed during the reaction process and 
we have distinguished between: 



1. Bonds between a free reactant and an active site (we 
will refer to the latter type of bond as single-bonds); 

2. Bonds between a linked reactant and an active site of 
another chain; 

3. Bonds between a reactant and two active sites of the 
same chain (which in the following we will call loops). 
We have analyzed the kinetics of bond formation for a 
system of size L = 100, with C r = 0.025 and n as = 5 
varying the chain concentration C and the probability 
Pb of bridge formation. The time is measured in Monte 
Carlo unit time St. 

Since C r -C C, the total number of bonds at the end 
of reaction is equal to twice the number of cross- linkers: 
N b = 2C r L 3 /8. In FigO the total number of bonds 
rib(t) is plotted as a function of the time together with 
the number of single-bonds n 8 (t) (bonds of type 1) and 
bridges n,b r (t) (i.e. bonds of type 2 or 3) in the case 
Pb = 0.01. The number of bonds has been normalized by 
Nb- As the reaction begins, single-bonds form rapidly, 
then bridges start to form and the degree of connectiv- 
ity between chains increases. The behavior of the total 
number of bonds versus time closely resembles the re- 
leased heat experimentally measured during the reaction 
reported in FigJ^l The velocity of the reaction is related 
to the probability pb of bridge formation. It governs the 
mean time of link formation between different chains and 
strongly influences the duration of the reaction process. 
In the inset of Fig[5] 7i&(i), n s (t) and rib r {t) are plot- 
ted as a function of the time for pb = 1, showing that 
both single-bonds and bridges form more rapidly as com- 
pared to the pb = 0.01 case. The total number of bonds 
nb(t) = n s (t) + 2ribr(t), and its time dependence can be 
written as n&(t) = iVj,(l — /(*))• In Fig[51 the function 
f(t) = 1 — nb(t)/Nb, representing the fraction of bonds 
that remain to form, is plotted in a semi logarithmic plot 
for the case pb = 0.01. The data reproduce extremely 
well the behavior observed experimentally (Fig|3]) and 
are well fitted by a sum of two exponentials 

f(t) =a 1 -e { - t ' T ^ +a 2 -e ( - t / T ^ (2) 

in agreement with the microcalorimetry measurements 
(EqfT]). From the fit we obtain t\ = 20 ± 2 and t 2 = 
166 ± 5, for pb — 0.01 and C — 0.3, providing t 2 /ti = 
8.3 ± 0.9. If the bridge probability pb varies, the mean 
time of bridge formation r 2 changes, and so does the ratio 
t%/t\. In the inset of Fig[5] f(t) is plotted for pb = 1. 
From the fit we obtain t\ = 34.4 ± 1.5, r 2 = 57.8 ± 1.5, 
and hence T 2 /n = 1.68 ± 0.08. 

In order to give a microscopic interpretation for these 
two characteristic times, we have directly computed the 
average time of formation of single bonds and bridges 
respectively, which cannot be done by microcalorimetry 
measurements. These two times are in agreement with 
the fitting parameters t\ and r 2 within error bars. Next, 
we have analyzed the mean square displacement of the 
reactant monomers when they are free (1) and when they 
have formed one single bond (2). Interestingly we have 
observed that the ratio between the corresponding diffu- 
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FIG. 5: Total number of bonds, number of bridges and num- 
ber of single-bonds normalized by the total number of possi- 
ble bonds Nb as a function of time for C = 0.3, CV = 0.025, 
n sa = 5 and pi, = 0.01. Inset: The same quantities for pb = 1. 
th indicates the time corresponding to half of the reaction. 



sion coefficients D\ /D2 is of the order of the ration T2 / T\ 
for all the analyzed concentrations of chains and cross- 
linkers. This result suggests that t\ and T2 are related to 
the characteristic times for diffusion of the free reactant 
and of the reactant linked to a gelatin chain, respectively. 

In agreement with this microscopic interpretation, our 
data (Fig[5] and O show that, for the concentrations C 
and C r explored, the single-bonds form more rapidly than 
bridges. This different velocity of formation is due, apart 
hompi,, to the different mobility of free cross-linkers with 
respect to linked ones, which are forced to move together 
to the chains to which they are permanently bonded. As 
a consequence, t\ < T2, i.e. the average time of formation 
of a single-bond is generally smaller than the average time 
of formation of bridges even for pb = 1 . 

It is interesting to notice that in simulations we can 
easily vary the bridge probability pb, to see how the fea- 
tures of the cross-linking reagent could possibly affect 
the kinetics of the bond formation. Remarkably, we find 
that this effect can be crucial. Indeed, as pb governs the 
bridges formation, it influences the average time T2: As 
Pb decreases, the reaction slows down and T2 increases. 
We have systematically analyzed the behavior of t^/ti, 
as a function of bridge probability pb- In FigJT] we plot 
the obtained data: When the bridge probability pb in- 
creases, the average time of bridge formation decreases, 
and so does the ratio ti/ti. For pb > 0.3 we find that 
the ratio decreases more slowly apparently tending to a 



FIG. 6: Function f(t) = 1 — rib(i)/Nb as a function of time, 
where Nb is the total number of possible bonds, for C r = 
0.025, n as = 5 and pt = 0.01. The full lines are the fitting 
curves ~ exp(—t/ri) with ri = 20, and ~ exp(—t/r2) with 
T2 = 166. Inset: f(t) with pb = 1. The full lines are the 
fitting curves ~ exp(—t/ri) with t\ = 34.4 and ~ exp(—t/T2) 
with T2 = 57.8. 

plateau. In this regime, the bridge probability pb does 
not influence the kinetics of bond formation, which is 
completely governed by the diffusion of the monomers 
and by their concentration. The value of the ratio be- 
tween the two characteristic times obtained in the exper- 
iments at pH = 6.7 (see Fig|3|) corresponds to a reactant 
with 0.005 < pb ^5 0.01 in our simulations. The exper- 
imental findings show that the solution pH also affects 
the kinetics of bond formation, i.e. the ratio t™/t™ de- 
creases as the pH decreases. This is in agreement with 
the chemistry of reaction, where the non protoned form of 
the amine is reactive. In fact, one could expect different 
regimes depending on the chains concentration. Indeed, 
if decreasing of the number of active sites will in general 
correspond to an increase of t\ , the effect of this variation 
on T2 is likely to strongly depend on chain concentration. 
At high concentration of chains, if the number of active 
sites per chain decreases below a certain level, we expect 
that these will be surrounded by many other sites of the 
same chain which are not active. As a consequence, due 
to excluded volume effects, they will hardly be reached 
by partially linked cross-linkers, i.e. T2 will increase as 
the number of active sites decreases and so will the ratio 
T-xjr\. On the other hand, T2 strongly depends on the 
chain mobility and on the formation of loops. Therefore 
decreasing the pH at low concentrations, leads to an in- 
crease of T\ that may be balanced by a not so dramatic 
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FIG. 7: The ratio tz/ti between the average time of forma- 
tion of bridges and the average time of formation of single- 
bonds as a function of bridge probability pb for different con- 
centrations C of chains, for CV = 0.025 and n as = 5. 



FIG. 8: The time of reaction end tf (dots), time of half 
reaction th (squares) and their ratio (triangles) as a function 
of chain concentration C for CV = 0.0025, n as = 5 and Pb = 
0.01. 



increase of T2, due to the eventual formation of loops. As 
a consequence one can observe a net decrease of the ratio 
t%It\. We would like to stress that, although the inves- 
tigation of the role of pH could be in principle done with 
this model by varying n as , one should use long enough 
chains to be able to detect the different concentration 
regimes. 

Conversely, when the number active sites increases up 
to n as — 10, in the analyzed concentration range the 
average time T\ does not change appreciably, while Ti 
decreases due to the formation of loops. For bridge prob- 
ability pb sufficiently high (p& > 0.8), the mean time of 
bridge formation becomes less than or equal to the mean 
time of formation of single-bonds: in these cases, the 
number of bonds Ub(t) may be fitted by a single expo- 
nential with the characteristic time T\. 

To complete our study of the kinetics of bond for- 
mation, we have measured the time th of formation of 
half of the total bonds, and the duration of the reac- 
tion tt, i.e. the average time needed to form all possi- 
ble bonds. The data are presented in Fig[5] and show 
that the ratio tf/th decreases as the concentration in- 
creases, tending to a plateau value for C > 0.3. This 
behavior is in agreement with the experimental findings: 
For Cbvsm = 0.3% g/cm 3 and C ge i ranging from 3% 
to 6%, the ratio tf/th decreases from tf/th ~ 27 to 
tf/th ~ 15. For Cg e i > 6% the ratio remains constant 
with the concentration. Comparing the experimental re- 
sults with simulations, we conclude that the regime of 



concentration C ~ 0.3 for the lattice model corresponds 
to C ge i ~ 6% of the experiments This correspon- 
dence is coherent with the behavior shown in Fig.s I5I6T71 
and support the suggested interpretation. It is worth to 
notice that t f /th 3> 1 for the explored range of parame- 
ters: this is a consequence of the fact that in the first half 
of the reaction most of formed bonds are single-bonds. In 
the second half, bridges between active sites are estab- 
lished, requiring a longer time. Finally we have analyzed 
the number of loops which are formed during the cross- 
linking reaction and may play a crucial role in the me- 
chanical response of the gel. Loops are not detectable in 
experimental measurements, but can be easily monitored 
in simulations. Previous numerical simulations of poly- 
merization process in hexamethylene diisocynate-based 
polyurethane [TtJ indicate that the number of loops has 
different roles in the various concentration regimes. In- 
deed the loss of elasticity due to loops may be outweighed 
by the increase of topological entanglements, depending 
on the concentration. 

At the end of the reaction we have counted the number 
of loops, normalized by the maximum number of bridges 
Nb/2, and investigated its behavior as a function of C. 
Data plotted in FigfS] refer to chain with n as = 5 and 
n as = 10. Our results indicate that the number of loops 
decreases as the chain density C increases. In the range 
of concentration explored the number of loops decreases 
following a power law behavior ~ C~ l characterized by 
an exponent I — 0.75 ± 0.05. The behavior appears in- 
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FIG. 9: Number of loops normalized by the maximum num- 
ber of bridges Nb/2, as a function of gelatin concentration 
C for p b = 0.002, 0.01, 0.1, 1 for n as = 5 and for Pb = 1 with 
rias = 10, for Cr = 0.025. The full line is the power law fitting 
curve ~ C~ l with I = 0.75. 

dependent of the bridge probability pb'- This result con- 
firms that the bridge probability only influences the ki- 
netics of bond formation but does not strongly affect the 
morphology of the system. By opportunely tuning pb, 
i.e. changing the reactant, the velocity of the reaction 
may be adjusted and the formation of single-bonds and 
bridges may be tuned in time, but the final geometrical 
properties of the structure should not be modified. On 
the other hand, the connectedness of the system may be 
influenced by the number of active sites per chain n as . 
In particular, when n as increases, the number of loops 
formed in the system increases and hence the degree of 
connectedness of the system decreases (FigJU). Moreover, 



in the limit of very diluted solutions of chains where all 
sites are active (n as = 10), loops represent about 80% 
of the total number of bonds and, as a consequence, the 
viscoelastic properties can be sensibly modified [ItJ ■ 



V. CONCLUSION 

In conclusion, our experimental results show that two 
different timescales affect the kinetics of bond forma- 
tion in our cross-linked gelatin solution. The numeri- 
cal data reproduce well the experimental ones and clar- 
ify the mechanisms involved in bond formation. Our 
study shows that the two time scales detected in exper- 
iments correspond to the average time of forming sin- 
gle bonds reactant-chains and bridges chains-chains via 
cross-linkers. These two times are related to the charac- 
teristic times of diffusion of free reactants and reactants 
which have already formed one bond. Their ratio controls 
the kinetics of the bond formation: Varying the concen- 
tration, the cross-linker reactivity and the pH strongly 
affect this ratio and therefore the kinetics of the gela- 
tion process. Our findings also show that the probability 
p b to form a bridge between two active sites allows to 
finely tune the kinetics of the reaction via the ratio of 
the two characteristic times. A variation of pb in our in- 
terpretation corresponds to a variation of the free energy 
barrier to be overcome in order to form the bond, or to 
different orientations of bonds vectors; hence to vary pb 
corresponds to change the reactant agent in the gelatin 
solution. Moreover, our data indicate that the number of 
loops formed between two active sites of the same chain, 
which has an important effect of the viscoelastic proper- 
ties of the system, increases when the pH of the solution 
increases. This model represents a useful tool to inves- 
tigate rheological behavior of gelatin solutions and the 
relation between the kinetics and gel structures. 
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